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Abstract. We propose a model for radio supernovae (RSN) 
based on the synchrotron emission from relativistic electrons 
which are diffusively accelerated at the expanding supernova 
shock. This model was originally developed for application to 



the optically thin emission observed from SN1987A. Here we 
generalise it by including the effects of free-free absorption 
I from both an external screen and from material internal to the 
| t source, and by relaxing the restriction to an azimuthal B-field. 
(^J We find a good fit to the entire set of radio data for the best 
observed highly-luminous RSN - SN1986J - with a reduced x 
of 3.85. Applying the new model to SN1988Z, another intrin- 
sically bright RSN, also yields a good fit (Xied ~ 2) but this 
is less significant, because of the limited data on this distant 
\£) (z = 0.02) source. These fits suggest that the shock expands 
£f) at constant speed, that the magnetic field within the source 
decreases with time according to t~ 2 , and that the compres- 
C ~ ) sion ratio of the shock front is close to the value expected of a 
r- H strong shock in an ideal gas of adiabatic index 5/3 - indicating 
a relatively low value of the cosmic ray pressure compared with 
ON SN1987A. In the case of SN1986J we derive an explosion date 
in August/September 1982, a magnetic field at the position of 
Q ^he shock 1000 days after explosion of B PS 4nT and a spatial 
I diffusion coefficient of the electrons of k fa 4 x 10 19 m 2 s -1 , 
O four orders of magnitude greater than the Bohm value. In ad- 
^__» dition, we obtain the optical depths to external and internal 
^ absorption, and derive an estimate of the mass-loss rate. 

Key words: acceleration of particles - shock waves - circum- 
stellar matter - supernovae: SN1986J, SN1988Z - radio con- 
tinuum: stars 



1. Introduction 

All radio supernovae (RSN) detected so far show a pattern 
typical for a source which is initially optically thick and later 
becomes optically thin: the emission appears first at high fre- 
quencies and extends later to lower ones. At fixed frequency 
the flux rises, reaches a peak and then fades somewhat more 
slowly than it rose. Once the radio emission is well past its 
peak, the flux density decreases with increasing frequency and 
is well-fitted by a power law of spectral index between 0.5 and 
1.2 - a strong indication that it is synchrotron radiation from 
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relativistic electrons. This behaviour indicates that the radio 
emission originates behind the outwardly propagating shock 
front. At early epochs the ionised component of the stellar 
wind of the progenitor upstream of the shock front is suffi- 
ciently dense for free-free absorption by thermal electrons to 
prevent the escape of the synchrotron radiation. As the shock 
propagates out through the overlying wind material the effec- 
tive optical depth decreases allowing the radiation to escape 
(Chevalier 1982; Weiler et al. 1986). 

However, the observed flux results from a convolution of 
the absorption process with the process responsible for the 
intrinsic emission, which itself depends on the way in which 
electrons are accelerated and on the evolution of the magnetic 
field within the source. Chevalier (1982) suggested that the 
intrinsic radio emission from RSN could be understood by as- 
suming that a fixed fraction (about 2%) of the flux of energy 
through the shock front is divided equally between relativistic 
electrons and magnetic field. Approximating the emission re- 
gion as a thin shell behind the shock, and including the effect 
of free-free absorption in a stellar wind with density falling 
off with the square of the distance from the progenitor star, 
Chevalier predicted that the radio flux density S(v, t) as a func- 
tion of frequency v and time t should be described by a five 
parameter model. Weiler et al. (1986) proposed a more gen- 
eral empirical model with six free parameters, of which the 
Chevalier (1982) model is a special case. The data from sev- 
eral Type lb and Type II RSN (most notably SN1979C and 
SN1980K) were shown by Weiler et al. (1986) to be well-fitted 
by both the general model and Chevalier's model. However, 
these models failed to fit one of the brightest and best-observed 
RSN: SN1986J. 

First discovered as a radio source, SN1986J was subse- 
quently found to have prediscovery radio and optical detec- 
tions. The explosion date remains uncertain; optical spectra 
indicate that it was a "peculiar" Type II supernova, and it is 
thought to have involved the explosion in 1982-83 of an ex- 
tremely massive progenitor with perhaps M ~ 20 — 60 Mq 
(Rupen et al. 1987; Weiler & Sramek 1988). The radio emis- 
sion from SN1986J rose much more slowly than is compatible 
with the rather rapid drop in absorption implied in the stan- 
dard model. This led to the suggestion that there might be an 
additional contribution to the free-free absorption due to ther- 
mal material which is mixed with the nonthermal synchrotron- 
emitting material (Weiler et al. 1990 - henceforth "WPS"), an 
effect which is easily incorporated for a slab-like geometry giv- 
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ing a significantly better fit to the data from early epochs. 

In previous work (Ball & Kirk 1992 - henceforth "BK") 
we developed a model for the unique second phase of radio 
emission from SN1987A (Staveley-Smith et al. 1992, 1993; Ball 
et al. 1995) based on the diffusive acceleration of electrons at 
the expanding supernova shock. It has since emerged that the 
emission from SN1987A requires for its understanding a de- 
tailed treatment of the environment of the progenitor and the 
reaction of accelerated particles on the hydrodynamics of the 
explosion (Duffy et al. 1995). In this Letter we suggest that in 
the case of SN1986J such complications can be neglected, and 
show that it is possible to fit a simple form of the model (includ- 
ing a combination of internal and external free-free absorption 
as proposed by WPS) to the radio data. This has the advan- 
tage that the free parameters are simply related to physical 
properties of the explosion, enabling estimation of quantities 
such as the compression ratio of the shock front, the magnetic 
field and the effective diffusion coefficient of the electrons. 

In Sect. 2, we describe the basic properties of the diffusive 
shock acceleration model. We identify the free parameters of 
the model and discuss their relation to the physical param- 
eters describing the shock and the environment into which 
it expands. In Sect. 3 we present the results of using a % 2 - 
minimisation procedure to fit the model to the long term evo- 
lution of the radio emission from SN1986J. In Sect. 4 we discuss 
this result and its implications, and note that fits to the obser- 
vations of SN1988Z, another peculiar Type II supernova, lend 
further support to the diffusive acceleration model. 
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where To = 3/(3 + 2a) and n = {1 + ^/[roit-l)]}- 1 , with 
Ri given by the solution of the transcendental equation 

v [R 1+T0 (i-i e )] (i+3n>/3 [3i? 1 /(2r a)]- 4/3 = 

f exp[3_Ri/(r a4cc)] • (3) 

[Note that the integral in Eq. (2) can be expressed in terms of 
incomplete Beta functions (BK).] This model for the intrinsic 
flux density has six free parameters all of which can be di- 
rectly identified with physical properties of the source and/or 
acceleration mechanism. They are: (i) the overall flux scaling 
parameter K\ ; (ii) the spectral index a which gives the shock 
compression ratio p c = (3 + 2a) / (2a); (iii) the epoch of ex- 
plosion t e ; (iv) the index n giving the time dependence of the 
magnetic field strength B oc (t — t e )~ n ; (v) the acceleration 
time iacc, which is related directly to the electron diffusion 
coefficient k and the compression ratio, provided the shock 
speed v B is known: k = t^ cc v 2 (p c — l)/[3p c (p c + 1)]; (vi) the 
frequency scaling parameter i>o, which is related to the mag- 
netic field strength at the shock front and the momentum po 
with which electrons are injected: vo = 10 3 "cii (po / m e c) 2 Biooo , 
where -Biooo is the magnetic field at the shock front at age 
1000 days, and m = 1.3 x 10 10 HzTesla" 1 . 



2. The model 

The model presented by BK assumes diffusive particle accel- 
eration at a shock expanding at constant speed so that the 
shock radius r s is proportional to the time after explosion, 
r s oc (t — t e ) where t = t/(lday) is the dimensionless time vari- 
able and i e is the time of the explosion. This is the behaviour 
expected in the similarity solutions found by Chevalier (1982) 
if the density profile of the ejecta is very steep. It is also consis- 
tent with almost all of the published fits to the radio emission 
of other RSN (Weiler et al. 1986). In the case of SN1987A, 
the magnetic field in the wind of the progenitor was assumed 
to be predominantly azimuthal. This leads to a time depen- 
dence of the magnetic field at the shock front B(r s ) of the 
form B(r s ) oc (i — f e ) _1 - If, on the other hand, the magnetic 
field encountered by the shock is predominantly radial, then 
B(r s ) oc (i — i e )~ 2 ■ Here we allow for each of these situations 
by setting B(r s ) oc (t — t e )~ n , where n is a free parameter. 
The model also assumes that the rate of injection of electrons 
into the acceleration mechanism is a constant: given that the 
shock speed does not vary, this corresponds to the injection of 
a fixed fraction of the electrons from a wind of density varying 
as the inverse square of the distance r from the progenitor. The 
resulting intrinsic (unabsorbed) flux density can be written as 
a function of the dimensionless frequency v = i//(5GHz) and 
time in the form 

F(v,t) = Ki v~ a (i — i e ) 1 ~ n ( 1 + a ^ ip(v,i,a,i e ,n,i^ cc ,v ) (1) 

where ij)(y, t, a, t e , n, i a cc, vo) is a well-behaved function, which 
is constant at large t, and which can be found from a straight- 
forward generalization of Eq. (8) of BK: 

ip(i>, i, a, i e , n, i^ cc , v ) = 
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Fig. 1. The intrinsic emission from the diffusive acceleration model 
at two frequencies differing by a factor of five — 1 GHz and 5 GHz 
say — for n = 1 and n = 2. 

Fig. 1 shows the intrinsic flux density as given by the diffu- 
sive acceleration model at two frequencies differing by a factor 
of five, for the cases n = 1 and n = 2. The plots clearly show the 
initial increase in the intrinsic flux density at a given frequency, 
the fact that the intrinsic emission looks "optically thin" at all 
times (i.e. the intrinsic flux density always decreases with in- 
creasing frequency), and that the eventual decay is a power-law 
in time which is steeper when n is larger. 
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The BK model for the second phase of radio emission from 
SN1987A does not include absorption effects because the on- 
set of the emission was optically thin, indicating that by then 
the shock had expanded sufficiently for the absorption to be 
negligible. This is not the case for the emission observed from 
other RSN so that the effects of absorption must be included. 
We treat free-free absorption in the same way as WPS, so that 
the radio flux density is given by 

S(v,t) = A e (T)Ai(T')F(v,i) (4) 

where F(y,t) is the intrinsic (unabsorbed) emission, 

A e (r)=e~ r with r = K 2 v~ 2 1 (i - i e )~ 3 (5) 

describes the external absorption by the wind, and 

Ai{r') = (f - e- T ')/r' with r = K z v~ 2 1 (t - i e )" 5 (6) 

describes the 'internal' absorption by material which is mixed 
with (at least in a statistical sense) the nonthermal radiating 
electrons. These formulae apply for a shock moving at constant 
speed. (Note that Eqs. (5) and (6) imply in the notation of 
WPS S = -3 and 8' = -5.) 

External absorption is related to the density of the environ- 
ment and, hence, to the mass-loss rate of the progenitor. The 
parameter if 2 is determined by the density of the progenitor's 
wind, which is assumed proportional to r~ 2 , and from Weiler 
et al. (1986), we find 

M = lO^/^Vo^Tf^Moyr- 1 , (7) 

where ww x lOkms -1 is the speed of the progenitor's wind, 
T4 x 10 4 K its temperature (after explosion) and v s = « s ,4 x 
lCkms" 1 . 

Internal absorption can be interpreted as due to either a 
homogeneous mixture of emitting and absorbing plasma, or to 
a statistical superposition of emitting and absorbing clumps 
along the line of sight. The parameter K3 is determined by 
the total number of thermal electrons mixed with the radiat- 
ing material, which is assumed constant. The relevant phys- 
ical parameter is the effective optical depth t[ 000 {5 GHz) at 
1000 days after explosion, measured at a frequency of 5 GHz: 
'"i'ooo(5 GHz) = 10 _15 ii' 3 - 

In total, the model we propose has eight free parameters: 
six for the intrinsic flux density and two more, if 2 and K3, 
for external and internal absorption. At early times the evo- 
lution of the flux density and the frequency spectrum are de- 
termined by a combination of the increasing intrinsic emis- 
sion and the decreasing absorption. At large times, when ab- 
sorption is negligible, the diffusive model gives a flux density 
S(0,i) oc v~ a (i - i e y- n(1 + a) so that for a given value of the 
spectral index a the flux decline at late times is faster when 
n is larger. (Note also that the asymptotic behaviour of the 
diffusive model is equivalent to that of the WPS model if in 
the notation of WPS one identifies fi = 1 — ra(l + a).) 

3. Fits to the observations 

We have used the ^-minimization technique to find the best fit 
of the diffusive acceleration model with both external and in- 
ternal absorption, as given by Eqs. (1-6), to the observed emis- 
sion from RSN. Using four-frequency observations of the ra- 
dio emission from SN1986J extending to 1991 June 11 (WPS; 
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Fig. 2. Best fit of the diffusive acceleration model (curves) to the 
observed radio emission from SN1986J up to 1991 June 11. Crosses 
indicate measurements at 1.5 GHz; circles 5 GHz; squares 15 GHz 
and triangles 22 GHz. 



Weiler, van Dyk & Sramek, private communication) we obtain 
an acceptable fit with \ 2 P er degree of freedom Xied = 3.85, 
where Xied = x I '(id — i p ) with Wd the number of data points 
and ra p the number of free model parameters. The parameters 
of this best fit model are given in Table f, and Fig. 2 shows 
the resulting flux densities together with the data. 

The value of Xied f° r the best fit diffusive model to the 
SNI986J data is comparable to that obtained for model fits 
of the Chevalier mini-shell model by WPS (to observations 
extending to 1988 December 28). The model parameters that 
are a direct analogue of those of the model fits by WPS are a, 
t e , if 2 and K3. The value of a that we obtain (a = 0.57) is 
somewhat smaller than that found by WPS (a = 0.67). This 
difference arises at least partly because the extended data set 
now available contains more information on the later optically 
thin stages of the emission. The value of t e quoted in Table f is 
expressed relative to the best fit explosion time found by WPS, 
and is consistent with their fit. The values of if 2 and K3 that 
we obtain are both considerably larger than the values given 
by the best fit model quoted by WPS, but differ somewhat in 
their definition, since WPS find a decelerating shock front. 

A fit of the diffusive acceleration model to the luminous ra- 
dio supernova SNf 988Z (see van Dyk et al. 1993 for a detailed 
discussion of the properties of this source), also gives an accept- 
able Xied ~ 2. The best-fit model implies that in SNI988Z, as 
in SNf986J, the magnetic field at the shock falls off with the 
inverse square of the time since explosion (B oc (i — i e )~ 2 ), 
and that the shock front is almost unmodified (p c fa 4). How- 
ever, the relatively small number of observations of this dis- 
tant source, and their necessarily large associated uncertain- 
ties, mean that the results of model fits to SNf988Z are less 
significant than those to SNf986J. Nevertheless, these sources 
are both peculiar Type II supernovae, which are qualitatively 
similar, and it is encouraging that the fits of the diffusive ac- 
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Table 1. Parameters for the best fit of the diffusive accelera- 
tion model to the observed radio emission from SNf986J. No 
uncertainties are quoted for K\, K2 and K3 because they are 
composites of the free parameters which were minimised. 



Parameter Physical interpretation 



Ki 


f.8xf0 lu 






f.f xfO 7 


M = 3.4 x fO-^io^Tf^Moyr" 1 




2.5xf0 15 


r i'ooo (5 GHz) = 2.5 (internal absorption) 


a 


0.58±0.0f 


p c = 3.6I ± 0.02 


te 


-2f±26 


t e = 1982 Aug. 24 


n 


2.0f±0.06 


B cx t~ 2 


t&cc 


332±25 


k = 4.2 x 10 19 ?> 2 4 m 2 s" 1 


fo/10 13 

Xred 


4.8±2.2 


-B1000 (po/m e c) 2 = 3.5 ± 1.7nT 


3.85 





celeration model imply that the physical source parameters are 
also similar. 

4. Discussion and Conclusions 

The salient points of the best-fit diffusive model to the radio- 
luminous SN1986J and SN1988Z are (i) a magnetic field de- 
pendence B oc t~ 2 and (ii) a compression ratio not too differ- 
ent from the value 4 expected at a shock front which is not 
modified by accelerated particle pressure. A magnetic field de- 
pendence of this type could arise if the field in the progenitor's 
wind were close to radial - perhaps as a result of a very small 
rotation rate of the star, or of turbulent processes within the 
wind itself. It is interesting to note that such a rapid decrease 
in the field at the position of the shock front is more difficult 
to achieve if the shock suffers appreciable deceleration, since 
in this case even a radial field would weaken more slowly than 
t~ . The value of the compression ratio indicates that the shock 
front has been slightly modified by the pressure of accelerated 
particles, but that this has not been as effective as in the case of 
SN1987A (Duffy et al. 1995). Such an effect may be connected 
with the sharp fall-off of the magnetic field. The pressure in 
relativistic particles is dominated by particles of TeV energy 
or higher, and a weak magnetic field will dramatically increase 
the time needed to accelerate particles to such energies. The 
magnetic field at the shock position 1000 days after explosion 
is found to be about 3 nT, significantly stronger than the am- 
bient interstellar field. This is consistent with our picture, in 
which the field is related to the progenitor. Extrapolating back 
to the stellar surface implies a prohibitively high surface field, 
but the model does not require the t~ 2 behaviour to apply 
earlier than about 500 days after explosion. The effective dif- 
fusion coefficient we derive is equivalent to a mean free path of 
an electron of 10 12 m, some four orders of magnitude greater 
than the gyro-radius of a radio emitting electron. This finding 
is similar to that for SN1987A (BK) although, in the present 
case, the magnetic field geometry does not require particles to 
propagate normal to the direction of the field. Our value of the 
mass-loss rate of the progenitor is somewhat larger than found 
by WPS, and the optical depth we find for free-free absorption 
by material mixed with the emitting source is slightly higher 
at 1000 days. However, as noted by WPS, the errors, especially 
in the estimate of the mass-loss rate, are large. 



The diffusive acceleration model also provides adequate 
fits to the radio emission from other less luminous but well- 
observed RSNe, namely the Type lb RSN SN1983N, and the 
"normal" Type II RSNe SN1979C and SN1980K. In contrast 
to SN1986J and SN1988Z, the fits to these other RSN imply 
a magnetic field dependence B oc t -1 which is consistent with 
the modelling of the second phase of emission from SN1987A 
(BK). This points to a significant difference in the progenitors 
of different supernova types. There are also indications that 
the emission from the recent RSN SN1993J requires a more 
detailed approach to the modelling, and in particular that the 
density of the circumstellar wind into which the shock is ex- 
panding decreases in this object more slowly than r~ 2 (van 
Dyk et al. 1994; Suzuki & Nomoto 1995). Thus, differences be- 
tween the models for the increasing number of well-observed 
RSN may well provide useful insights into the properties of 
these sources, and help provide a physical basis for the various 
supernova classifications. 

The diffusive model has more free parameters (a total of 
eight) than the Chevalier mini-shell model and the more gen- 
eral empirical models previously used to interpret the observa- 
tions of RSN. However, the value of x P er degree of freedom 
provides an objective measure of the quality of the fit, inde- 
pendent of the number of free parameters, and the best fits of 
the diffusive acceleration model are comparable in quality to 
the fits of the other models to SN1986J (WPS), SN1988Z (van 
Dyk et al. 1993) and to less luminous RSN such as SN1979C 
(Weiler et al. 1986). Furthermore, the proposed model has the 
advantage that it involves a specific physical mechanism for 
accelerating electrons to the relativistic energies required to 
produce the observed radio emission. 
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